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1.  Introduction 


Reliable  wireless  communication  is  of  paramount  importance  for  many  important 
civilian  and  military  applications.  In  scenarios  where  the  propagation  environment 
is  challenging,  such  as  in  complex  terrain,  urban  canyons,  and  indoor  environ¬ 
ments,  achieving  a  reliable  communications  link  that  is  not  significantly  impeded 
by  the  presence  of  obstacles  and  scatterers  is  extremely  challenging.  In  recent  years, 
rapid  global  urbanization  and  the  ever-increasing  need  for  reliable  communications 
in  emergency  and  tactical  applications1,2  has  necessitated  reliable  communications 
paradigms  for  such  challenging  environments.  Networks  are  needed  to  support  hu¬ 
mans  and  autonomous  agents,  incorporating  communications,  sensing,  real-time 
positioning,  and  tracking.3,4 

State-of-the-art  communications  systems  that  operate  in  the  upper  very  high  fre¬ 
quency  (VHF)  through  microwave  bands  suffer  from  significant  attenuation  and 
small-scale  fading,  as  well  as  phase  distortion  caused  by  high  levels  of  multipath.  5~7 
In  a  complex  environment  at  these  frequencies,  without  line  of  sight  (LOS),  mul¬ 
tipath  propagation  is  the  dominant  signal  source  between  a  transmitter  and  a  re¬ 
ceiver,  whereas  the  direct  path  signal  is  often  extremely  weak  due  to  penetration 
loss.  Major  effort  has  been  devoted  over  many  decades  to  design  sophisticated  mul¬ 
tipath  mitigation  techniques,  incorporating  channel  coding  and  modulation  schemes 
to  overcome  the  challenges  introduced  by  the  channel.8^13  Such  schemes  add  sig¬ 
nificantly  to  transceiver  hardware  and  software  complexity,  still  may  only  provide 
marginal  performance  in  challenging  scenarios,  e.g.,  necessitating  networking  over 
many  hops  or  a  high  transmit  power. 

The  lower  VHF  band  has  significant  potential  to  support  low  complexity,  low  power, 
highly  reliable  communications.  This  primarily  stems  from  the  fact  that  in  the  lower 
VHF  scatterers  are  small  in  terms  of  wavelength.3,14  Consequently,  strong  penetra¬ 
tion  through  multiple  walls  and  buildings  can  be  achieved  at  relatively  low  power. 
Reflection,  scattering,  and  diffraction  phenomena  are  dramatically  reduced,  thereby 
greatly  minimizing  multipath  fading,  yielding  a  short-range  channel  that  is  LOS- 
like  in  terms  of  very  slight  phase  distortion  and  delay  spread.  This  liberates  the 
system  designer  at  low  VHF  from  the  typically  stringent  requirements  on  power, 
system  bandwidth,  and  complex  equalization  processing  needed  in  ultra-high  fre¬ 
quency  (UHF)  and  microwave  based  systems. 
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Of  course,  conventional  antennas  such  as  dipoles  operating  in  the  lower  VHF  band 
are  very  large,  limiting  application  of  such  systems  to  stationary  cases  or  large  plat¬ 
forms.  In  the  case  of  mobile  nodes  with  smaller  antennas,  the  transmit  power  must 
be  high  enough  to  compensate  for  the  poor  efficiency  of  the  antennas.  Thus  far, 
the  large  size  of  conventional  antennas  has  been  a  major  obstacle  in  terms  of  their 
large-scale  deployment  in  mobile  scenarios  within  highly  cluttered  environments. 
However,  due  to  recent  advances  in  antenna  miniaturization  techniques  and  the  de¬ 
velopment  of  palm-sized  lower  VHF  antennas  with  good  performance,1516  interest 
in  low  power,  low  data  rate  communications  in  this  band  is  increasing.  Thus,  there 
is  a  need  to  better  understand  the  channel  characteristics  in  this  band,  so  as  to  better 
inform  and  guide  the  design  of  systems  operating  at  the  lower  VHF  band. 

While  the  frequency  allocations  vary  in  different  regions  of  the  world,  in  North 
America,  there  is  dual  allocation  at  low  VHF,  consisting  of  primary  and  secondary 
allocations.  For  example,  the  primary  allocation  near  38  and  40  MHz  is  for  Federal 
use,  including  research  and  development  and  test  range  operations.  Other  applica¬ 
tions  and  users  include  the  industrial,  scientific,  and  medical  (ISM)  band,  the  US 
Army’s  Single  Channel  Ground  and  Airborne  Radio  System  (SINCGARS),  high¬ 
way  patrol  radios,  and  ocean  radars.1718 

The  large  majority  of  channel  characterization  in  complex  propagation  scenarios 
has  consisted  of  path  loss  measurements  in  support  of  long-range  commercial  wire¬ 
less  applications,  which  operate  in  the  upper  UHF  and  microwave  range. 19-24  These 
studies  typically  assume  a  powerful  base  station  with  an  elevated  antenna  serving 
mobile  units  located  near- ground, 1  as  is  common  in  cellular  phone  and  other  mobile 
communications.  In  contrast,  little  attention  has  been  paid  to  shorter-range  lower 
VHF  band  channel  characterization  for  near-ground  nodes  deployed  in  complex 
propagation  scenarios,  especially  for  indoor- to-indoor  and  outdoor-to-indoor  com¬ 
munications.  Urban  path-loss  models  in  the  upper  VHF  and  military  UHF  bands 
based  on  field  measurements  have  been  presented.1'25  Similar  path-loss  models  in 
the  lower  VHF  range  have  also  been  proposed.2'26-33  However,  classical  empirically 
derived  urban  propagation  models,  often  geared  toward  conventional  cellular  ap¬ 
plications,  in  which  a  base  station  is  situated  at  a  relatively  high  elevation,  are  not 
appropriate  for  near- ground  nodes  and  such  models  are  not  necessarily  valid  at  the 
lower  VHF  and  high  frequency  (HF)  bands.31'32  Here  we  primarily  focus  on  measur¬ 
ing  the  channel  transfer  function  and  phase  distortion  in  the  lower  VHF;  extensive 
path-loss  measurements  have  been  reported,33  including  experimental  comparison 
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of  path  loss  at  40  MHz  and  2.4  GHz  for  indoor  and  indoor/outdoor  cases. 

Simulation-based  propagation  models  have  also  been  proposed  to  understand  chan¬ 
nel  characteristics  of  multipath  scenarios  in  VHF/UHF  bands.  34  36  In  previous  stud¬ 
ies,37  4(1  investigations  of  short-range  HF/VHF  propagation  in  forested  environ¬ 
ments  have  been  presented.  There  has  also  been  some  effort  in  developing  systems 
operating  in  the  HF  band  for  near-field  communication  and  geolocation  applica¬ 
tions.41,42  We  note  that  the  majority  of  research  regarding  propagation  and  channel 
modeling  in  the  HF/VHF  band  is  concerned  with  long-range  ionospheric  propaga¬ 
tion.43,44 

In  this  report,  we  study  near-ground,  wireless  channel  modeling  in  the  lower  VHF 
band  for  indoor  and  indoor/outdoor  scenarios,  drawing  from  extensive  propagation 
measurements  and  simulations  using  a  finite-difference  time-domain  (FDTD)  based 
electromagnetic  model.  In  Section  2,  the  full-wave  electromagnetic  simulation  re¬ 
sults  are  described.  In  Section  3,  our  40-MHz  measurement  system  is  discussed, 
including  our  software-defined  radio  (SDR)  transmit  and  receive  tools,  as  well  as 
integration  onto  a  robotic  platform  for  ease  of  data  collection  over  large  areas.  In 
Sections  4  and  5,  we  present  the  measurement  scenarios,  experimental  results,  and 
analysis  based  on  our  measurements,  including  path  loss,  small-scale  fading,  phase 
distortion,  and  bit  error  rate  (BER)  measurements. 

2.  Multipath  Analysis  and  Full-Wave  Simulation 

In  this  section,  we  consider  the  multipath  channel  and  associated  channel  transfer 
function.  This  is  a  precursor  to  measuring  and  estimating  the  channel  transfer  func¬ 
tion  in  amplitude  and  phase.  An  indoor/outdoor  propagation  simulation  is  used  as 
an  initial  look  at  the  problem  (later  we  consider  a  physical  scene  that  is  very  similar 
to  the  simulation).  These  results  set  the  stage  for  our  experiments  in  the  sections 
that  follow,  including  statistical  tests  for  the  presence  or  absence  of  multipath,  test¬ 
ing  for  an  ideal  transfer  function,  and  communications  BER  performance  without 
equalization. 

Our  simulation  geometry  is  shown  in  Fig.  1 .  In  the  upper  UHF  to  microwave  bands, 
ray  tracing  approaches  are  often  used  to  simulate  wave  propagation.  Ray  tracing 
is  relatively  straightforward  and  computationally  tractable  under  the  assumption 
that  the  size  of  the  scatterers  in  the  environment  is  generally  much  larger  than  the 
wavelength.  However,  this  assumption  breaks  down  in  the  lower  VHF,  where  the 
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size  of  some  scatterers  becomes  comparable  to  the  wavelength.  For  this  reason, 
we  resort  to  a  technique  that  directly  solves  Maxwell’s  equations  to  characterize 
the  wireless  channel  in  the  lower  VHF  band.  In  particular,  we  use  EMCUBE,  a 
commercial  full-wave  solver  based  on  FDTD.45 


Fig.  1  The  FDTD  simulation  geometry  used  for  short-range  phase  distortion  analysis.  The 
walls  are  modelled  as  homogeneous  dielectric  slabs  and  metal  sheets  are  introduced  to  simulate 
the  propagation  characteristics  of  windows. 

The  simulation  scene  in  Fig.  1  consists  of  several  walls,  floor,  and  ceiling,  as  well  as 
doors  and  windows  (the  floor  and  ceiling  are  not  shown  in  the  figure  for  clarity).  The 
floor  is  modeled  as  a  half-space  homogeneous  dielectric  medium  with  a  dielectric 
constant  of  concrete  (er  =  4.5  +  j  0.0 11).  The  walls  are  modeled  as  dielectric  slabs 
with  a  dielectric  constant  of  brick  (er  =  4  +  jO.OOl).  To  make  this  indoor  scenario 
more  realistic,  rectangular  metal  sheets  ( 1 .2  m  by  1 .2  m)  are  included  to  simulate  the 
effect  of  windows  as  shown  in  the  figure.  The  dimensions  of  the  overall  scenario  are 
15  m  by  15  m  by  3.3  m.  The  geometry  is  kept  relatively  small  in  order  to  make  the 
simulation  more  computationally  tractable.  A  vertically  polarized  dipole  antenna 
positioned  outdoors  is  used  as  a  transmitter  (Tx),  emitting  a  Gaussian  envelop  pulse 
that  is  frequency  translated  to  40  MHz,  with  a  1-MHz  bandwidth.  Four  vertically 
polarized  receive  (Rx)  antennas  are  positioned  indoors  and  the  received  signal  from 
each  antenna  is  recorded.  Note  that  the  Rx  positions  are  symmetrically  positioned 
on  either  side  of  the  Tx  antenna  along  the  x  axis. 

In  addition  to  path  loss,  we  consider  the  phase  coherence  and  memory  in  the  chan¬ 
nel.  Denote  the  linear  time  invariant  channel  impulse  response  as 
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(1) 


N—l 

h(n)  =  Y  pk^kS(n  -  nk), 

k=0 


where  Pk  and  <pk  are  the  amplitude  and  phase  of  the  kth  component  of  the  signal. 
The  k  —  0  term  corresponds  to  the  first  arrival  in  time,  where  we  assume  no  <  n\  < 
■  ■  ■  <  un- i  model  the  delays  of  the  N  multipath  components.  Taking  the  Fourier 
transform,  we  can  express  the  transfer  function  as 


N- 1 

H(f)  =  P0ej*0e-j2nfn 0  +  Y  Pkej*ke-j2nfn\ 

k= 1 


(2) 


where  Pk  and  <f>k  are  the  amplitude  and  phase  of  the  kth  component  of  the  Fourier 
transform  of  h(n).  The  transfer  function  phase  can  be  written  as 


lH(f)  =  tan  1 

Generally,  in  the  presence  of  multipath,  the  transfer  function  phase  in  Eq.  3  is  a 
nonlinear  function  of  frequency.  If  the  amplitude  of  secondary  multipath  compo¬ 
nents  is  small,  the  first  term  in  Eq.  2  will  be  the  dominant  component.  In  that  case, 
all  terms  k  >  1  in  Eq.  2  can  be  ignored,  and  only  the  leading  term  in  the  numerator 
and  denominator  of  Eq.  3  will  be  significant,  resulting  in  a  phase  that  is  linear  in 
frequency  given  by  0  —  2n fn0.  In  Fig.  2,  the  phase  of  the  received  signal  at  each 
receive  antenna  position  is  plotted  against  frequency,  centered  around  40  MHz.  The 
linear  phase  variation  (aside  from  phase  wrapping)  in  the  bandwidth  of  interest  is 
evidence  of  the  minimal  effect  of  multipath  in  this  scenario.  Similar  results  (not 
shown)  occur  in  the  transfer  function  amplitude,  which  is  flat  across  the  simulated 
bandwidth.  Although  there  is  attenuation  due  to  penetration,  there  is  little  to  no  ev¬ 
idence  of  multipath  propagation  or  fading.  For  example,  in  the  case  of  receiver  po¬ 
sition  R2,  there  are  4  layers  of  walls  between  the  transmit  antenna,  which  is  outside 
the  building  and  the  receiver  located  indoors.  The  received  signal  is  predominantly 


P  sin(0o  -  2nfnp)  +  J2k=i  pk  sin(^fc  ~  27r  fnk) 
P  cos(0o  -  27T fnQ)  +  J2k=i  pk  cos(0fc  -  2irfnk) 
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the  component  that  penetrates  through  the  walls  along  the  shortest  path  between 
the  Tx  and  Rx  antennas.  Reflections  and  multiple  scattering  between  walls  are  very 
small  due  to  the  small  size  of  the  walls  (note  that  the  thickness  of  a  typical  wall  is  a 
small  fraction  of  the  wavelength  at  40  MHz).  If  there  was  significant  multipath,  the 
phase  variation  would  be  nonlinear,  which  is  not  the  case,  as  can  be  seen  in  Fig.  2. 


Fig.  2  Phase  variation  of  the  received  signal  as  a  function  of  frequency  for  the  4  Rx  positions 
Ri  to  i?4  as  shown  in  Fig.  1,  computed  via  an  FDTD  full-wave  electromagnetic  simulation  of 
the  scene. 

3.  Measurement  System 

In  this  section,  we  describe  the  hardware,  peripherals,  and  software  that  we  use 
to  experimentally  collect  and  analyze  data.  Figure  3  illustrates  the  2  measurement 
systems  employed.  The  2  systems  are  essentially  identical,  except  that  in  Fig.  3b 
we  incorporate  a  reference  channel  by  simultaneously  transmitting  the  signal  via 
the  antenna  and  via  cable  to  the  receiver,  enabling  waveform  distortion  comparison, 
whereas  in  Fig.  3a  the  signal  is  only  received  via  an  antenna. 

3.1  System  Hardware 

3.1.1  Radio  and  Daughtercards 

We  use  the  USRP  model  N210  SDR  for  our  tests.46  The  USRP  N210  provides 
an  Ethernet  interface  that  allows  us  to  use  software  control  for  transmission  and 
reception.  The  USRP  is  comprised  of  2  parts:  a  motherboard  consisting  of  an  field- 
programmable  gate  array  (FPGA),  analog-to-digital  convertor  (ADC),  and  digital- 
to-analog  convertor  (DAC),  and  a  modular  daughtercard  whose  functions  vary  by 
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model,  but  may  include  filtering  and  up-  and  down-frequency  conversion.  The  USRP 
N210  has  a  100  MS/s  ADC  and  400  MS/s  DAC.  We  use  the  BasicTx  and  BasicRx 
daughtercards  at  the  transmitter  and  receiver,  respectively,  that  provide  access  to 
signals  in  the  1  to  250  MHz  range.  These  daughtercards  do  not  have  local  oscilla¬ 
tors  or  up/down  conversion  capability;  we  carry  out  these  operations  in  the  software. 
In  our  tests,  we  use  2  USRPs,  one  a  dedicated  transmitter  and  the  other  a  dedicated 
receiver. 


RF  cable 


Fig.  3  USRP  based,  1-way  (dual  channel)  measurement  system  diagrams  for  a)  pulse  and 
packet  tests,  and  b)  tone  tests.  The  pulse  and  packet  tests  are  used  to  characterize  the  channel 
phase  distortion  and  BER,  respectively.  For  the  tone  test,  a  reference  channel  is  introduced  by 
cable  connection,  allowing  direct  comparison  of  the  wireless  and  wired  signals. 


3.1.2  Antenna 

The  antenna  used  for  these  experiments  is  a  A/6  short  dipole  whose  center  fre¬ 
quency  is  tuned  to  40  MHz.  Because  short  dipoles  have  high  capacitive  reactance, 
a  pair  of  high-Q  loading  coils  are  inserted  on  both  arms  of  the  dipole  to  increase 
the  reactance  and  hence  cancel  out  the  capacitive  effect,  resulting  in  better  match¬ 
ing  with  the  transmission  line  feeding  the  antenna.  The  antenna  has  a  gain  of  -2 
dBi  and  a  fractional  bandwidth  of  0.075.  The  antenna  was  measured  using  the  out¬ 
door  elevated  antenna  measurement  technique  by  mounting  the  antennas  on  40-ft 
fiberglass  towers.  Simulated  Sll  and  radiation  pattern  of  the  antenna  are  shown  in 
Fig.  4.  While  in  this  report  we  focus  on  channel  measurement  and  characteriza¬ 
tion,  evolving  miniature  antenna  designs  operating  in  the  upper  HF  and  lower  VHF 
band,  with  dimensions  of  A/100  and  smaller,  have  significant  potential,  especially 
for  short-range,  low  power  applications.15 

3.1.3  Mobile  Robotic  Collection  Platform 

To  efficiently  conduct  large-scale  experiments,  we  integrated  our  USRP  receiver 
onto  a  robotic  platform,  enabling  data  collection  over  extended  indoor/outdoor  envi- 
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Frequency  [MHz] 

Fig.  4  Simulated  Sll  and  radiation  pattern  of  the  A/6  antenna  used  for  all  measurements.  The 
center  frequency  was  tuned  to  40  MHz  using  the  coils  on  the  antenna. 

ronments.  The  robot  has  the  capability  to  carry  out  a  simultaneous  localization  and 
mapping  (SLAM)  algorithm  to  provide  a  map  and  locate  itself  within  the  map.  For 
our  experiments,  the  robot  is  teleoperated  throughout  each  scenario.  By  integrating 
the  USRP  with  the  robot  control  system,  each  radio  measurement  is  automatically 
supplemented  with  time  and  position  information.  This  enabled  collection  of  large 
datasets  over  extended  space  and  time  (for  example,  see  Scenario  3  in  Section  4). 

Referring  to  Fig.  5,  we  use  the  iRobot  PackBot,  outfitted  with  a  Hokuyo  UTM- 
30LX  laser  scanner,  a  pan-tilt  unit,  global  positioning  system  (GPS),  an  onboard 
computer  with  Solid-state  drive  (SSD),  and  a  Microstrain  3DM-GX2  inertial  mea¬ 
surement  unit  (IMU);  all  of  these  instruments  along  with  robot  odometry  data  may 
be  input  to  the  SLAM  algorithm  for  localization.  The  overall  position  error  at¬ 
tributable  to  the  system  is  about  ±1.5  m.  To  this  baseline  setup,  we  integrated  a 
USRP  N210  SDR,  the  A/6  antenna,  and  a  BB-2590  military  battery  cradle  that 
supplies  power  to  the  SDR. 

A  Python  script  running  on  the  onboard  computer  interfaces  with  the  robot  operat¬ 
ing  system  (ROS),  where  the  robot  control  and  SLAM  algorithms  reside.46^9  This 
allows  us  to  teleoperate  the  robot  using  an  Xbox  360  wireless  controller.  When  the 
robot  stops,  the  Python  script  invokes  a  GNU  radio  function  that  captures  a  user- 
specified  interval  of  data  at  a  specified  center  frequency  and  sampling  rate.  The 
radio  transmitter  continuously  emits  the  channel  sounding  signal  of  interest  (tones 


8 


or  pulses,  depending  on  the  particular  experiment).  The  sampled  radio  data,  as  well 
as  GPS,  odometry,  IMU,  and  laser  scan  data  from  the  robot,  are  then  recorded  for 
later  post-processing. 


Laser 

Scanner 


BB-2590  Battery 
and  Cradle 


Packbot 


Antenna 

USRP 

Fig.  5  The  SDR  receiver  system  integrated  on  a  packbot  mobile  platform.  Note  the  antenna  is 
only  partially  visible. 

3.2  Data  Acquisition  Software 

We  use  signal  processing  software  from  the  open-source  GNU  Radio  project50  and 
a  custom  MATLAB-USRP-Library  (MUL)51  to  control  the  USRPs  and  conduct  our 
experiments.  GNU  Radio  enables  the  implementation  of  SDR  functions  including 
filters,  modulators,  demodulators,  channel  coders,  and  many  other  elements. 

For  our  tests,  we  suitably  modify  GNU  Radio’s  built-in  code  for  our  purposes;  this 
code  includes  pure  tone  generation  as  well  as  complex  baseband  sample  acquisition 
of  the  transmitted  waveform.  The  MUL  is  a  software-library  that  connects  MAT- 
LAB  to  the  USRP,  allowing  rapid  conversion  of  communications-oriented  MAT- 
LAB  scripts  into  actual  wireless  experiments.  We  implemented  a  packetized  com¬ 
munications  link  through  MUL  that  allows  for  the  generation  and  capture  of  pack¬ 
etized  data  that  we  use  for  the  BER  measurements  described  in  Section  4.2.2. 
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4.  Measurements  and  Processing 


4.1  Measurement  Scenarios 


Measurements  were  carried  out  in  3  different  environments,  2  indoor  (Scenarios 
1  and  2)  and  1  extended  indoor/outdoor  (Scenario  3).  Different  cases  are  listed  in 
Table  1,  and  experimental  parameters  are  shown  in  Table  2.  The  Tx  and  Rx  positions 
in  Table  2  correspond  to  those  marked  in  Figs.  6  through  8,  as  described  next. 

Table  1  Measurement  scenarios  indicating  transmit  locations  and  receiver  regions,  with  ref¬ 
erence  to  Figs.  6, 7,  and  8,  respectively.  Most  measurements  consisted  of  NLOS  propagation. 


Environment 

Tx  position 

Rx  locations 

Setup 

Case  la 

Scenario  1 

Pi 

P.1-P0 

NLOS 

Case  lb 

Scenario  1 

Pi 

P7  -  Pl3 

NLOS 

Case  lc 

Scenario  1 

Pi 

P3 

NLOS 

Case  2a 

Scenario  1 

P2 

P3-P5 

NLOS 

Case  2b 

Scenario  1 

Po 

Pi 

LOS 

Case  3a 

Scenario  2 

Pi 

P2-P7 

LOS 

Case  3b 

Scenario  2 

Pi 

P0-P7 

LOS 

Case  3c 

Scenario  2 

Pi 

P7  -  PU 

NLOS 

Case  3d 

Scenario  2 

Pi 

P8-P9 

NLOS 

Case  4 

Scenario  2 

Po 

Pl2  ~  Pl6 

NLOS 

Case  5a 

Scenario  3 

Po 

Ro 

NLOS 

Case  5b 

Scenario  3 

Pi 

f?4  -  R7 

NLOS 

Case  5c 

Scenario  3 

Po 

i?8,-Rg 

NLOS 

Case  5d 

Scenario  3 

Pi 

R  \  -  R3 

NLOS 

Case  5e 

Scenario  3 

Po 

Rs,R9 

NLOS 

Table  2  Parameters  used  for  measurements 


Parameter 

Setting 

Max.  Tx  power 

15  dBm  (32  mW) 

Tx,Rx  Ant.  Gain 

-6.5  dBi,  -6.5  dBi 

Tx  Ant.  height 

1  m(0.13A) 

Rx  Ant.  height 

0.8  m,  1  m  (0.1A) 

Tx,  Rx  pol. 

v,v 

Tx  Amp.  Gain 

40  dB 

Scenario  1  consists  of  indoor  measurements  in  a  multi- story,  house-like  test  facility. 
Figure  6  shows  a  schematic  of  the  first  floor  of  this  building,  which  has  the  same 
geometry  as  the  one  simulated  (see  Fig.  1  in  Section  2).  This  is  a  wooden  structure, 
including  the  floors,  walls,  and  ceilings.  As  indicated  in  Fig.  6,  steel  cylinders  and 
sheets  were  inserted  as  might  occur  in  a  home  environment. 


10 


P7  P8 
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Fig.  6  Measurement  Scenario  1:  2nd  floor  layout  of  wooden  building,  with  steel  sheets  and 
cylinders  introduced  as  shown 


Scenario  2  is  depicted  in  Fig.  7,  interior  to  an  office  and  laboratory  building  with 
extensive  metal  partitions,  reinforced  walls,  and  metal  objects.  Measurements  were 
taken  in  the  hallways,  adjacent  to  offices  and  large  interior  lab  bays.  Measurements 
in  Scenarios  1  and  2  were  taken  up  to  a  maximum  range  of  40  m,  in  LOS  and  NLOS 
conditions. 

For  Scenario  3,  the  test  facility  shown  in  Fig.  8  was  used,  with  a  test  diameter 
of  about  200  m,  consisting  of  multiple  buildings  constructed  of  concrete  and  cin- 
derblock,  with  fences  and  other  dielectric  and  metallic  objects  in  the  scene.  Thou¬ 
sands  of  measurements  were  collected  for  the  2  different  transmitter  positions  de¬ 
noted  PO  and  PI  in  the  figure,  teleoperating  the  mobile  robotic  collection  platform 
within  the  denoted  oval  collection  regions  and  along  the  dashed-line  collection  path 
shown. 
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13m 


i co 


Fig.  7  Measurement  Scenario  2:  An  L-shaped  corridor  in  the  2nd  floor  of  a  much  larger  office 
and  laboratory  building  with  substantial  metallic  partitions  and  walls 


Fig.  8  Measurement  Scenario  3:  overhead  view  of  the  test  facility.  P0  and  Pi  are  the  2  trans¬ 
mitter  positions  and  R0  to  R9  are  various  indoor  and  outdoor  regions  traversed  by  the  robot 
for  data  collection 
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4.2  Pulse  Tests 


In  this  section,  we  describe  our  pulse-based  test  for  the  presence  of  phase  distortion 
(caused  primarily  by  multipath  propagation).  These  experiments  are  based  on  peri¬ 
odically  transmitting  root-raised  cosine  (RRC)  pulses,  a  ubiquitous  waveform  used 
in  digital  communications.52 

4.2.1  Pulse  Test  Parameters 

Let  pin)  denote  the  discrete-time  RRC  pulse.  Figure  9  shows  pin)  and  its  corre¬ 
sponding  analog  waveform  pit)  after  passing  through  the  DAC.  Referring  to  the 
RRC  formula  in  Proakis,52  pit)  has  roll-off  factor  /3  =  0.9,  symbol  time  Ts  = 
3.2  x  10~5  s,  and  pulse  duration  1.28  x  10~4  s.  We  transmit  a  periodic  pulse  train 
with  8.72  x  10-4  s  off-time  between  pulses  during  which  null  samples  are  fed  to  the 
DAC.  Thus,  the  total  time  between  2  consecutive  pulse  starts  is  T  =  1.28  x  10~4  + 
8.72  x  HP4  =  HP3  s.  We  transmit  in  both  the  in-phase  (I)  and  quadrature  (Q) 
channels,  denoting  the  complex  baseband  equivalent  as  pc{t)  =  pit)  +  jp(t).  The 
RRC  pulse  has  a  10  dB  bandwidth  of  approximately  30  kHz,  comfortably  within 
the  bandwidth  of  our  antenna. 


Fig.  9  The  root-raised  cosine  waveform  used  in  our  pulse  tests.  Blue  stems  indicate  digitally 
generated  samples  at  the  transmitter  and  the  red  curve  the  actual  waveform  output  by  the 
DAC  for  transmission.  Pulse  duration  is  1.28  x  lO-4  s. 


In  our  tests,  the  transmitter  location  is  fixed  and  the  pulse  train  is  continuously 
emitted.  The  receiver,  mounted  on  a  robot,  is  navigated  to  a  given  point  and  stopped. 
Data  are  then  collected  for  |  s,  after  which  collection  is  terminated  until  the  receiver 
navigates  to  the  next  collection  point.  For  each  unique  receiver  position,  125  pulses 
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are  collected,  corresponding  to  the  transmitted  complex  baseband  signal: 


125 

Pr(t)  =  ^2 pc(t-kT ), 
k= 1 


(4) 


where  T  =  10~3  and  the  receiver  samples  the  passband  signal  at  100  MSPS,  deci¬ 
mates  the  sampled  signal  to  25  MSPS,  and  then  digitally  mixes  the  signal  to  base¬ 
band  using  a  local  estimate  of  the  carrier  frequency.  The  resulting  digital  signal  at 
complex  baseband  is  denoted  by 


r(n)  =  pr(n)  *  h(n)  +  b(n),  (5) 

where  h(n)  is  the  channel  impulse  response,  and  bin)  is  assumed  to  be  an  additive 
white  Gaussian  noise  (AWGN)  random  process.  In  the  analysis  that  follows,  we 
assume  that  the  channel  is  constant  for  each  |  s  collection  time,  i.e.,  all  125  pulses 
see  an  identical  channel. 

The  received  complex  baseband  pulse  train  r(n)  is  passed  through  a  matched  filter 
(MF)  whose  template  is  the  RRC  pulse  shape.  Denote  by  pm  the  output  of  match 
filtering  pc(n)  against  p(n);  note  that  pAf(n)  is  a  complex  valued  signal  with  the 
nth  sample  having  the  same  value  on  both  the  I  and  Q  channels.  We  then  estimate 
timing  and  compute  the  average  of  the  MF  individual  pulse  responses,  given  by 


PrH 


1 

125 


125 

5^PAf,*(rc), 

i=  1 


(6) 


where  pM,i(n )  denotes  the  MF  response  to  the  ith  received  pulse.  Given  the  RRC 
pulse  shape,  PM,iip)  is  ideally  a  raised  cosine  shape  if  there  were  no  noise  or  chan¬ 
nel  distortion.  The  additive  noise  effect  is  minimized  by  averaging  the  individual 
MF  responses. 

To  test  for  deviation  from  an  ideal  channel  without  multipath,  we  compute  the  ob¬ 
served  channel  transfer  function  given  by 
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(7) 


Hr  = 


Hpr) 

HpmY 


where  T  denotes  the  discrete  Fourier  transform  (DFT)  computed  via  an  fast  Fourier 
transform  (FFT)  algorithm.  Here  P(pr)  is  the  DFT  of  the  measured  averaged  re¬ 
ceived  pulse  after  matched  filtering,  and  T{pm )  is  the  DFT  of  the  ideal  matched 
filter  response  (a  raised  cosine  pulse).  Without  multipath,  the  LTI  channel  consists 
of  scaling  and  a  single  delay  only,  as  discussed  in  Section  2.  We  normalize  the  ini¬ 
tial  delay  to  zero,  and  the  maximum  transfer  function  amplitude  to  unity,  so  that 
ideally  \Hr\  =  1  with  phase  equal  to  zero. 

Multipath  propagation  induces  frequency-dependent  deviations  from  the  ideal  trans¬ 
fer  function  in  both  magnitude  and  phase,  so  we  consider  2  measures  of  channel 
distortion.  Let  S  denote  the  selected  set  of  frequencies  under  test,  taken  from  the 
DFT.  We  define  a  mean-square  error  (MSE)  of  the  difference  between  the  measured 
and  ideal  unity  magnitude  response  given  by 


M  =  T^[|ffK(/)|-1]2.  (8) 

11  fes 

The  MSE  of  the  difference  between  the  measured  channel  phase  response  and  its 
ideal  zero  value  is  given  by 


p  =  (9) 
11  feS 

These  2  measures,  sample  averages  computed  across  the  DFT  frequency  bin  set  S, 
capture  the  discrepancy  of  the  measured  channel  HR  from  an  idealized  all-pass, 
zero-phase  filter.  Note  that  for  each  fixed  transmit-receive  location,  there  is  a  sin¬ 
gle  channel  realization  Hr,  and  hence  one  value  for  M  and  one  for  P  for  each 
experimental  receiver  position. 

Recall  our  decimated  sampling  rate  is  25  MSPS  with  a  pulse  bandwidth  of  30  kHz, 
so  we  are  considerably  above  the  Nyquist  rate.  We  use  a  1-sided  DFT  frequency 
grid  of  0  to  50  kHz,  with  grid  spacing  A /  «  6  Hz.  The  set  S  consists  of  those 
frequency  values  on  the  grid  where  the  spectrum  of  the  ideal  received  raised  cosine 
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pulse,  |  J-’(pm)  |,  exceeds  0.01,  where  the  maximum  magnitude  has  been  normalized 
to  1.  Small  values  for  P  provide  evidence  of  minimal  phase  distortion,  and  small 
values  for  M  are  indicative  of  a  constant  attenuation  with  respect  to  frequencies 
probed  by  the  pulse.  As  described  in  the  next  section,  in  Scenario  3,  we  employ 
the  robotic  platform  to  repeat  the  above  procedure  for  hundreds  of  measurement 
locations  (keeping  transmitter  location  fixed  while  varying  receiver  location). 

4.2.2  Wider  Bandwidth  Pulse  Tests 

The  transmit  pulse  bandwidth  can  easily  be  increased  or  decreased  in  the  SDR 
transmitter  by  changing  the  sampling  rate  when  generating  the  RRC  waveform.  In 
addition  to  the  above  case,  we  also  generated  an  RRC  pulse  stream  with  symbol 
duration  of  T$  =  1.935  x  10~6  with  single  pulse  duration  of  7.74  x  10-6  s.  In 
this  case  we  transmit  10,000  pulses  per  second,  so  that  the  total  time  between  2 
consecutive  pulse  starts  is  T  =  10-4  s.  The  RRC  pulse  has  3-dB  bandwidth  of 
approximately  500  kHz  and  a  10-dB  bandwidth  of  800  kHz.  The  1-sided  DFT  grid 
over  which  we  evaluate  M  and  P  for  this  case  ranges  from  0  to  1  MHz.  Receiver 
processing  for  this  case  then  follows  as  described  above  for  the  narrowband  pulse 
tests. 

4.3  Tone  Tests 

A  second  phase  distortion  analysis  technique  is  also  pursued,  which  is  based  on 
tones  transmitted  simultaneously  over  the  air  and  via  a  cabled  reference  channel, 
shown  in  Fig.  3b.  These  tests  were  conducted  by  moving  the  receiver  to  each  new 
Rx  location.  Tone  frequencies  are  selected  within  [39.5,  40.5]  MHz.  The  transmit 
signal  is  divided  between  the  antenna  and  an  SMA  cable  (reference)  using  a  power 
divider.  Likewise  at  the  receiver,  complex  baseband  samples  of  the  2  separate  sig¬ 
nals  are  acquired  simultaneously  from  the  2  channels  by  using  the  2  ports  on  the 
BasicRx  card,  which  are  controlled  by  the  same  clock.  For  a  given  frequency,  the 
instantaneous  received  phase  is  computed  for  both  the  antenna  signal  and  the  refer¬ 
ence.  The  difference  in  instantaneous  phase,  averaged  over  all  samples  collected  at 
a  given  frequency  /,  yields  our  estimate  of  the  phase  response  at  /. 

4.4  Bit  Error  Rate  (BER)  Measurements 

In  addition  to  the  above  channel  sounding  tests,  BER  measurements  were  carried 
out  at  a  center  frequency  of  40  MHz,  where  the  receiver  estimates  timing  and  carrier 
offset,  but  does  not  equalize  for  any  delay-spread  in  the  channel.  At  the  transmitter, 
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packets  of  length  1000  symbols  are  generated  with  random  data,  quadrature  phase 
shift  keying  (QPSK)-modulated  with  Gray  coding,  and  per  symbol  RRC  filtered. 
We  transmit  62,  000  QPSK  symbols  per  second,  corresponding  to  a  data  rate  of  120 
kbps.  A  Barker  code  prefix  is  also  incorporated  into  each  packet  to  aid  synchroniza¬ 
tion.53 

At  the  receiver,  the  complex  baseband  signal  is  sampled  at  1  MSPS.  After  applying 
an  RRC  pulse  matched  filter,  symbol  timing  is  recovered  using  the  early-late  gate 
method.54  Next,  cross-correlation  against  the  Barker  code  prefix  is  used  to  estimate 
and  correct  the  frequency  offset.  Finally,  the  Barker  code  is  also  used  to  estimate 
and  correct  for  phase  offset  and  attenuation  via  a  complex  scalar  multiplication.  We 
emphasize  that  equalization  consists  of  a  single  scalar  complex  multiplication,  and 
no  attempt  is  made  to  estimate  or  compensate  for  any  presumed  delay  spread  in 
the  channel.  The  QPSK  BER  is  estimated  over  108  bits,  and  the  observed  signal-to- 
noise  ratio  (SNR)  recorded. 

5.  Experimental  Results 

5.1  Path  Loss 

We  first  report  on  the  path  loss  observed  in  Scenario  3.  We  use  the  pulse  data  de¬ 
scribed  in  Section  4,  for  Case  5a  in  Table  1  (outdoor,  NLOS).  For  each  fixed  Rx 
location,  raw  I  and  Q  samples  are  collected  spanning  125  pulses,  which  are  used 
to  compute  an  uncalibrated  pulse  power  (Pruncal).  It  is  uncalibrated  because  the  way 
the  USRP  processes  data,  the  I  and  Q  samples  do  not  map  linearly  to  the  actual 
receive  power;  USRP  gain  calibration  is  discussed  below.  Simultaneously,  we  also 
store  GPS  location  information.  The  GPS  data  are  post-processed  and  enhanced  for 
accuracy.  Such  enhancement  is  accomplished  by  leveraging  the  auxiliary  sensors  of 
the  robot,  including  IMU,  laser  scanner,  and  odometry,  to  localize  the  platform  with 
reference  to  a  fixed  Universal  Transverse  Mercator  (UTM)  origin.  We  then  correct 
the  processed  GPS  coordinates  by  introducing  a  translation  offset  using  measured 
ground  truth  and  Google  Maps  as  a  reference. 

To  calculate  the  path  loss  at  a  fixed  location,  the  out-of-band  noise  is  reduced  with 
a  narrowband  bandpass  filter.  Next,  Pruncal  is  computed  by  finding  the  peak  value 
of  |Pc(n)|2  (the  “peak  power”)  for  each  pulse  in  the  pulse  train  received  at  that 
location,  and  then  averaging  across  all  these  peak  power  values.  After  the  peak 
power  values  are  calculated,  the  channel  path  loss  (in  dB)  is  computed  using 
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(10) 


where  Pt  and  Pruncal  are  the  transmit  and  receive  powers  in  dBm,  respectively.  Here, 
Gt  and  Gr  are  the  gains  of  the  Tx  and  Rx  antennas,  and  Gsys  is  the  sum  of  all  other 
system  gains  and  losses  including  amplifier  gain,  filter  loss,  and  cable  losses.  Gsys 
also  includes  a  calibration  factor  for  the  USRP  to  map  /  and  Q  values  to  true  power 
levels.  This  calibration  factor  was  computed  by  directly  connecting  the  USRP  to  a 
signal  generator  via  an  radio  frequency  (RF)  cable  and  performing  a  power  sweep 
to  map  the  quantity  captured  by  (Pruncal)  to  actual  signal  power  in  dBm. 

Next,  the  corrected  GPS  data  are  used  to  find  the  distance  between  the  receiver  and 
transmitter  at  all  measurement  points.  The  range  data  are  binned  into  intervals  of 
size  0.45A,  with  the  path-loss  data  mapped  to  its  corresponding  bin,  and  the  mean 
and  standard  deviation  of  all  path  loss  values  falling  in  a  given  bin  are  computed. 
The  binning  size  is  experimentally  chosen  to  be  similar  to  the  position  error  of 
our  system.  The  overall  position  error  attributable  to  our  measurement  system  is 
within  ±1.5  m,  which  comes  from  errors  in  the  GPS  data  and  SLAM  processing.  It 
should  be  noted  that  the  linear  path  loss  trend  does  not  change  when  the  bin  size  is 
reduced,  although  the  path-loss  curve  will  not  be  as  smooth  and  has  a  slight  increase 
in  standard  deviation  for  the  linear  fit. 

The  path-loss  mean  and  standard  deviation  versus  range  is  shown  in  Fig.  10,  which 
is  based  on  measurements  in  Scenario  3  (case  5a)  and  consists  of  a  total  of  152 
collection  points.  The  path  loss  for  this  case  has  a  standard  deviation  of  2  dB  from 
the  least  squares  fit.  The  calculated  path-loss  exponent  is  6.3.  This  path-loss  expo¬ 
nent  value  is  slightly  higher  than  typical  indoor/outdoor  path-loss  exponent  values 
reported  in  the  literature  for  the  upper  VHF,  UHF,  and  microwave  range  (roughly 
4  to  6. 5). 55,56  In  the  lower  VHF  with  low  antenna  heights  (<  0.5A),  propagation 
is  primarily  supported  by  surface  waves,  resulting  in  the  higher  path-loss  exponent 
value  because  surface  waves  decay  faster  in  general.36'57  Note,  however,  that  while 
the  path-loss  exponent  is  relatively  high,  the  overall  attenuation  is  much  less  than 
expected  at  higher  frequencies  due  to  lower  penetration  loss  through  the  walls.33 
It  should  also  be  noted  that  other  measurements  in  Scenario  3  (with  different  Tx 
locations)  have  resulted  in  path-loss  exponent  values  ranging  from  4  to  5.  The  path- 
loss  variation  depends  on  the  number  of  obstacles  and  buildings  in  the  path.  Our 
reported  measurements  in  Scenario  3  are  a  challenging  case  with  a  significant  num- 
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ber  of  structures  in  the  typical  LOS  path. 


Fig.  10  Mean  (dashed  line)  and  standard  deviation  (solid  vertical  bars)  of  measured  path  loss 
in  Scenario  3  (Fig.  8  Case  5a  in  Table  1).  The  solid  line  is  the  least  squares  lit.  The  path-loss 
statistic  shown  is  computed  using  a  total  of  152  measurement  points  with  over  10,000  measured 
pulses,  where  each  pulse  consists  of  3,200  samples. 

A  more  extensive  measurement-based  path-loss  characterization  for  VHF  indoor/outdoor 
channels,  as  well  as  comparison  against  measurements  at  2.4  GHz,  is  presented  in 
a  previous  paper.33  Although  the  measurements  were  carried  out  in  different  in¬ 
door  and  outdoor  scenarios,  the  path-loss  results  presented33  are  comparable  to  the 
path-loss  results  discussed  in  this  report. 

5.2  Phase  Distortion  and  Instantaneous  Phase 

Multipath  or  other  propagation  effects  can  generally  introduce  phase  distortion.  To 
look  more  closely  at  this,  we  examine  the  instantaneous  phase  for  2  representative 
pulse  cases.  In  Fig.  11  we  plot  instantaneous  phase  against  time,  for  2  complex 
baseband  pulses  from  Scenario  3  case  5b,  receiver  location  Re  (see  Table  1  and  Fig. 

8).  The  instantaneous  phase  plotted  in  Fig.  11  consists  of  the  middle  1250  samples 
of  the  received  pulse,  out  of  3200  samples  per  pulse.  Receiver  processing  consisted 
of  bandpass  filtering  to  reduce  out-of-band  noise,  passband  sampling,  frequency 
offset  estimation,  and  mixing  to  complex  baseband.  No  matched  filter  processing 
or  averaging  was  used.  The  plots  represent  arbitrarily  chosen  60th  and  120th  pulses 
of  the  125  collected  at  this  particular  Rx  position,  and  the  SNR  was  measured  to  be 
close  to  20  dB  for  each  pulse.  In  fact,  the  entire  ensemble  of  125  pulses  all  show  the 
same  instantaneous  phase  as  in  Fig.  11.  Note  that  quadrature  results  in  a  constant 
noise-free  pulse  phase  of  7t/4,  which  has  been  removed  for  plotting.  Also,  carrier 
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frequency  offset  error  will  result  in  a  linear  phase  change  with  time;  we  estimate 
and  remove  this  linear  trend  before  plotting. 


Instantaneous  phase  at  baseband 


Fig.  11  Instantaneous  phase  vs.  time  for  2  example  complex  baseband  received  pulses  after 
frequency  offset  correction,  for  a  receiver  position  selected  from  Case  5b  (see  Fig.  8  and  Table 
1).  Little  phase  distortion  is  evident  for  high  SNR  cases  such  as  these. 

After  removing  the  7t/4  offset  and  linear  trend  due  to  carrier  offset,  plots  of  instan¬ 
taneous  phase  display  remarkably  little  variation  from  zero  for  these  and  all  other 
high  SNR  cases  examined,  demonstrating  very  little  observed  phase  distortion.  The 
slight  variation  in  the  phase  plots  is  well  within  what  can  be  expected  due  to  the 
inherent  limitations  of  our  measurement  system. 

5.3  Statistical  Analysis  of  the  Measured  Channel  Transfer  Function 

Next  we  report  on  the  measured  statistics  of  the  channel  transfer  function  as  de¬ 
scribed  in  Section  4,  using  the  error  measures  M  and  P  from  Eqs.  8  and  9,  respec¬ 
tively  (Figs.  12-14).  To  establish  baseline  distributions  for  M  and  P  for  each  fixed 
SNR,  we  simulate  a  channel  that  introduces  only  attenuation  and  AWGN  (an  ideal 
all-pass  zero-phase  transfer  function).  We  perform  the  identical  transmit  and  re¬ 
ceive  data  generation  and  processing  that  we  do  in  our  actual  data  collection  based 
on  processing  a  pulse  train,  as  described  in  Section  4.  For  each  SNR  value,  we  con¬ 
ducted  500  Monte  Carlo  trials  over  realizations  of  the  additive  noise.  We  used  the 
resulting  empirical  distributions  of  M  and  P  to  find  a  lower  1  —  a  confidence  inter¬ 
val  for  each,  i.e.,  we  find  threshold  values  tm  and  rP  such  that  Pr(M  >  tm)  =  ol 
and  Pr(P  >  t>)  =  a.  Note  that  these  threshold  values  are  a  function  of  the  test 
SNR.  For  our  tests,  we  use  a  =  0.05  to  yield  a  95%  confidence  level. 
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Fig.  12  Mean  square  deviation  of  the  received  signal  phase  compared  with  an  ideal  channel 
transfer  function  for  cases  5a  and  5b  (see  Fig.  8  and  Table  1).  P  measures  deviation  from  the 
ideal  zero  channel  transfer  function  phase.  Little  statistical  deviation  is  observed  at  the  95% 
level. 
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Fig.  13  Mean  square  deviation  of  the  received  signal  amplitude  compared  with  an  channel 
transfer  function  for  Cases  5a  and  5b  (see  Fig.  8  and  Table  1).  M  measures  deviation  from  the 
ideal  unity  gain  transfer  function  amplitude.  Only  a  small  deviation  is  observed. 


o 

Scenario 

3 (Case  5b) 

□ 

Scenario 

3  (Case  5a) 

+ 

Baseline 

statistic (no  multipath) 

o 

□ 

□ 

□ 

q 

n  °  h 

0  D 

n  D  o 

□  F 
8 


°Dd°°  ° 


O 

>°  □  n 


<*><*> 

o 

O 


<> 

"8 


o 

0°6o  T 
O  op 


I 


20 

SNR  [dBl 


21 


3.5 

3 

2.5 

"  2 
-a 

£ 

1.5 

1  - 

0.5- 


x  10 


Simulated  multipath  _ 
(two  path  model) 


I 


o 

> 

X 


o 

0.05,10 

* 

0.05,20 

• 

0.05,30 

□ 

0.1,10 

❖ 

0.1,20 

V 

0.1,30 

> 

0.15,10 

< 

0.15,20 

+ 

0.15,30 

X 

Measured 

+ 

Baseline 

statistic 

10  15  20 

SNR  [dB] 


25 


30 


Fig.  14  Mean  square  deviation  of  the  measured  phase  for  case  3d  (see  Fig.  8  and  Table  1), 
along  with  various  realizations  of  theoretical  multipath  channels  based  on  a  2-path  channel 
model.  The  first  and  second  numbers  in  the  legends  of  the  theoretical  multipath  points  are  the 
normalized  amplitude  and  delay  of  the  secondary  (multipath)  component  with  respect  to  the 
dominant  component.  Note  that  the  dominant  component  is  normalized  to  have  energy  1  and 
delay  0.  A  95%  confidence  interval  for  the  distribution  of  the  phase  statistic  under  an  ideal  flat 
fading  channel  with  AW GN  is  shown  (blue  bar).  Note  that  most  multipath  channels  fall  outside 
the  blue  bar;  this  shows  our  test  has  good  statistical  power  in  being  able  to  detect  multipath. 
Also  note  that  actually  obtained  experimental  values  (green  x)  falls  within  or  very  near  the 
blue  bar. 
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Figures  12  and  13  depict  experimental  results  for  measurements  carried  out  in  Sce¬ 
nario  3  and  Scenario  1 .  These  results  consist  of  1 1  and  97  different  Rx  positions 
in  Scenario  1  and  Scenario  3,  respectively.  The  data  collected  at  each  measurement 
point  consists  of  250  pulses  in  the  case  of  Scenario  1  and  125  pulses  in  the  case  of 
Scenario  3  resulting  in  a  total  of  14875  measured  pulses.  For  each  SNR,  the  inter¬ 
vals  (0,  tm)  and  (0,  tp)  are  also  shown  as  vertical  blue  bars.  Note  that  in  Fig.  13, 
these  bars  have  a  slight  upward  bias  from  their  ideal  lower  end  point  of  0.  This  is 
attributable  to  a  small  departure  from  unity  gain  due  to  roll-off  in  the  magnitude 
response  of  the  lowpass  filter  used  in  processing  the  pulses.  Statistical  evidence  of 
deviation  from  the  ideal  transfer  function  exists  at  level  a  if,  for  a  particular  com¬ 
puted  value,  M  >  tm,  or  P  >  tp. 

The  phase  statistic  P  in  Fig.  12  shows  that  all  the  Rx  positions  considered  resulted 
in  a  channel  transfer  function  whose  phase  deviation  from  an  ideal  delay  and  atten¬ 
uation  channel  is  well  within  statistical  limits  caused  by  additive  noise.  We  cannot 
reject  the  ideal  transfer  function  model  at  the  a  =  0.01  level.  The  amplitude  statistic 
shown  in  Fig.  13  shows  similar  results,  with  most  measurements  within  the  statis¬ 
tical  limits,  although  there  is  some  deviation.  However,  this  statistical  deviation  is 
small,  and  appears  to  be  independent  of  the  phase  results.  A  more  refined  analysis 
could  consider  the  joint  distribution  of  M  and  P,  although  only  a  small  variation  in 
attenuation  is  observed,  and  the  transfer  function  phase  seems  more  pertinent  to  the 
measurement  of  phase  distortion. 

The  phase  and  amplitude  statistics  are  also  computed  for  the  wideband  pulse  tests 
described  in  Section  4.  Figure  14  shows  the  phase  statistic  calculated  based  on 
measurements  carried  out  in  Scenario  2  where  a  few  LOS  and  NLOS  points  were 
measured.  Similar  to  the  previous  case,  the  phase  deviation  in  2  of  the  3  cases  falls 
within  the  95%  confidence  intervals  of  the  null  hypothesis;  the  third  is  just  slightly 
outside  it.  In  this  figure,  we  also  show  the  theoretical  phase  MSE  for  various  2- 
component  multipath  channels,  i.e.,  channels  consisting  of  a  dominant  component 
and  a  single  multipath  component.  We  show  channels  whose  secondary  component 
has  energy  which  is  a  fraction  a  of  the  dominant  signal,  and  delay  of  7  meters 
relative  to  the  dominant  signal,  for  various  values  of  a  and  7  (as  shown  in  the 
figure  legend).  For  example,  the  point  a  =  0.05,  7  =  10  represents  the  phase 
deviation  statistic  value  we  would  obtain  for  a  2-component  multipath  channel, 
whose  secondary  component  is  0.05  as  strong  as  the  primary  one,  and  which  travels 
an  excess  distance  (in  free  space)  of  10  m.  Note  that  for  the  higher  SNR  cases, 
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most  multipath  channels  give  phase  MSE  values  that  fall  outside  the  null  hypothesis 
interval;  this  indicates  that  the  test  has  strong  statistical  power.  Also,  note  that  in 
general,  as  the  secondary  component’s  relative  energy  and/or  delay  increases,  the 
value  of  the  statistic  grows.  This  shows  that  our  statistic  P  has  the  desirable  property 
that  it  grows  larger  as  the  multipath  becomes  worse  in  amplitude  and  delay.  Finally, 
note  that  the  various  multipath  channel  statistics  are  point  estimates  computed  in 
the  absence  of  AWGN.  Had  we  also  introduced  AWGN,  the  points  would  become 
intervals  and  would  shift  strictly  up,  as  AWGN  increases  phase  MSE.  We  have  not 
shown  these  latter  intervals,  so  as  not  to  clutter  the  figure. 

The  amplitude  statistic  shows  qualitatively  similar  results;  however,  the  Tx  and  Rx 
antennas  exhibit  small  amplitude  variations  (i.e.,  are  not  constant  gain)  over  the 
pulse  bandwidth,  which  affects  our  magnitude  MSE  statistic.  Note  that  this  is  more 
of  an  issue  for  the  wider  bandwidth  pulses  than  the  narrow  bandwidth  pulses,  as  the 
larger  bandwidth  results  in  a  larger  overall  deviation  from  constant  gain.  The  net 
effect  of  this  system  non-ideality  is  to  add  noise  to  our  magnitude  MSE  statistic. 
Rather  than  carefully  compensate  for  the  effect  of  this  antenna  gain  variation  on 
our  magnitude  statistic,  we  have  focused  on  the  phase  statistic. 

5.4  Tone  Tests  for  Transfer  Function  Phase  Measurement 

Next  we  present  results  for  tone  tests  as  described  in  Section  4.2.2,  where  we  mea¬ 
sure  the  phase  difference  between  the  received  tone  via  the  antenna  and  via  the 
cabled  reference  channel.  Experimental  results  are  plotted  in  Figs.  15,  16,  and  17, 
with  the  specific  cases  listed  in  the  captions  of  each  figure.  The  curves  show  the  re¬ 
sults  for  NLOS  cases  collected  in  Scenario  1,  shown  in  Fig.  6.  At  each  measurement 
location,  we  collect  2  x  106  (I,Q)  samples  from  the  2  channels  for  each  frequency. 
With  the  sampling  rate  of  2  Ms/s,  this  corresponds  to  a  1-s  collection  time.  The 
phase  difference  is  calculated  for  each  sample,  and  the  mean  and  standard  deviation 
of  the  phase  difference  is  plotted. 

The  experimental  results  show  a  strong  linear  phase  relationship  as  a  function  of 
frequency  for  the  entire  range  of  test  tones  in  [39.5,  40.5]  MHz.  Some  deviation 
is  expected  due  to  varying  SNR  and  other  experimental  effects,  but  the  measured 
deviation  is  small  as  shown  in  the  figures  (note  the  vertical  scale  is  in  radians).  This 
approach  is  a  direct  way  of  measuring  the  presence  or  absence  of  phase  distortion, 
since  phase  deviation  is  measured  with  respect  to  the  low  distortion  cable  reference 
channel.  The  observed  low-deviation  linear  phase  as  a  function  of  frequency  con- 
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firms  that  the  channel  largely  consists  of  a  single  dominant  propagation  path  for  all 
cases  tested. 


Case  2a:  NLoS  Case 
4. 5  | - t - 
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Fig.  15  Tone  tests:  Mean  and  standard  deviation  of  the  measured  phase  difference  between 
the  wireless  and  reference  channels  for  Case  2a,  at  Rx  point  P3  (see  Fig.  6  and  Table  1) 
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Case  2a:  NLoS  Case 
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Fig.  16  Tone  tests:  Mean  and  standard  deviation  of  the  measured  phase  difference  between 
the  wireless  and  reference  channels  for  Case  2a,  at  Rx  point  P4  (see  Fig.  6  and  Table  1) 


We  note  that  phase  and  amplitude  distortion  will  occur  as  the  frequency  increases, 
an  effect  that  will  become  evident  within  the  VHF  band.  In  a  previous  publica¬ 
tion,  the  authors  investigated  the  propagation  characteristics  at  HF/lower  VHF  as 
compared  to  100  MHz  and  higher,  using  full-wave  simulations.  3  For  a  given  indoor 
building  scenario,  the  effect  of  multipath  and  channel  induced  distortion  is  clearly 
evident  at  100  MHz,  but  is  very  minimal  at  lower  frequencies  such  as  20  or  40 
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MHz.  This  is  easily  seen  in  the  signal  (both  amplitude  and  phase)  coverage  inside 
the  building.3  A  single  tone  test,  when  carried  out  at  several  different  locations  in  a 
relatively  confined  environment,  will  show  evidence  of  multipath  in  both  amplitude 
and  phase,  if  such  multipath  is  present. 


Case  2a:  NLoS  Case 
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Fig.  17  Tone  tests:  Mean  and  standard  deviation  of  the  measured  phase  difference  between 
the  wireless  and  reference  channels  for  Case  2a  at  Rx  point  P5  (see  Fig.  6  and  Table  1) 


5.5  BER  Measurement  Results 

Finally,  we  report  on  BER  tests,  as  described  in  Section  4.2.2.  These  tests  were 
carried  out  in  Scenario  2,  with  the  Tx  location  fixed  and  the  Rx  position  varied 
as  shown  in  Fig.  7,  yielding  both  LOS  and  NLOS  cases.  Figs.  18  and  19  show  the 
measured  BER  vs.  SNR  curves  for  the  LOS  and  NLOS  scenarios,  respectively.  Also 
plotted  is  a  theoretical  QPSK  BER  for  an  AWGN  channel.  The  measured  SNR  was 
estimated  for  each  Rx  position  (generally  declining  as  range  was  increased).  We  see 
excellent  agreement  between  theoretical  and  measured  BER,  with  some  small  devi¬ 
ations  attributable  to  finite  size  effects  arising  from  the  number  of  bits  collected.  As 
noted  in  Section  4.2.2,  the  receiver  processing  only  corrected  for  carrier  frequency 
offset  and  phase  correction,  and  did  not  attempt  to  correct  for  any  delay  spread.  Un¬ 
like  at  higher  frequencies  such  as  higher  VHF,  UHF,  and  microwave,  the  observed 
lower  VHF  channel  in  a  dense  indoor  scene  with  large  metal  obstacles  behaved 
as  an  ideal  AWGN  channel,  such  that  ideal  BER  results  were  obtained  with  very 
simple  receiver  processing. 
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Fig.  18  Theoretical  and  measured  BER  vs.  SNR  curves  for  LOS  indoor  channel  Scenario  2 
(Case  3b,  See  Table  1  and  Fig.  7).  Simple  receiver  processing  achieves  the  theoretical  perfor¬ 
mance  limit  predicted  for  an  AWGN  channel. 
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Fig.  19  Theoretical  and  measured  BER  vs.  SNR  curves  for  NLOS  indoor  channel  Scenario  2 
(Case  3d,  See  Table  1  and  Fig.  7).  The  theoretical  performance  limit  is  attained  in  this  NLOS 
case. 
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6.  Conclusion 


In  this  report,  we  have  characterized  the  lower  VHF  band  wireless  channel  around 
40  MHz  in  highly  cluttered  environments  using  an  extensive  measurement  cam¬ 
paign  including  statistical  analysis  of  thousands  of  transmissions.  Our  objective 
was  to  study  the  viability  of  low  power,  short-range,  reliable  communications  that 
are  minimally  affected  by  the  presence  of  clutter,  and  establish  the  degree  to  which 
phase  distortion  is  present  or  absent  as  is  relevant  to  geolocation  in  these  environ¬ 
ments. 

FDTD  electromagnetic  simulations  of  an  indoor/outdoor  scene  showed  little  phase 
distortion  or  multipath  propagation  that  is  so  characteristic  of  upper  VHF  and  above. 
This  helped  motivate  our  study  of  phase  distortion,  channel  transfer  function,  and 
communications  BER.  Tone  and  pulse-based  measurements  were  conducted  in  3 
different  scenarios,  primarily  in  NLOS  cases,  over  ranges  of  10  to  200  m.  Using 
SDRs  for  both  transmit  and  receive,  and  integrating  the  receiver  onto  a  robotic  plat¬ 
form,  gave  us  experimental  flexibility  to  collect  tens  of  thousands  of  measurements. 

Shorter-range  measurements  carried  out  inside  2  different  buildings  (Scenarios  1 
and  2),  as  well  as  medium-range  measurements  performed  in  a  complex  indoor/outdoor 
propagation  environment  (Scenario  3)  were  used  to  quantitatively  study  path  loss, 
phase  distortion,  and  channel  transfer  function.  It  was  shown  that  the  channel  phase 
distortion  is  minimal  both  in  LOS  and  NLOS  scenarios,  even  when  several  layers 
of  obstacles  are  present  between  the  Tx  and  Rx  nodes.  The  channel  phase  distortion 
was  analyzed  using  2  independent  techniques,  including  tone  tests  that  incorporated 
a  cable  reference  channel  for  precise  phase  comparisons. 

Pulse  tests  were  used  to  estimate  the  empirical  channel  transfer  function  in  [39.5,  40.5] 
MHz,  and  we  found  no  statistical  evidence  against  a  simple  scalar  attenuation  and 
AWGN  channel  model.  We  performed  tone  tests  by  which  the  ideally  expected  lin¬ 
ear  variation  of  phase  as  a  function  of  frequency  was  confirmed.  We  also  showed 
that  BER  measurements  carried  out  in  the  same  environments  show  very  good 
agreement  with  the  theoretically  predicted  AWGN  channel  BER  curve,  for  both 
LOS  and  NLOS  cases. 

Thus,  our  central  finding  is  that  for  low  power,  short-range  low  data  rate  commu¬ 
nications  in  cluttered  environments,  the  lower  VHL  channel  can  be  modeled  as  an 
AWGN  channel  with  scalar  attenuation.  The  simplicity  of  the  channel,  along  with 
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recent  advances  in  the  design  of  extremely  miniaturized  lower  VHF  antennas,  can 
be  exploited  in  a  large  variety  of  signal  processing  and  communications  applications 
including  geolocation  in  GPS -denied  environments  and  voice  communications  in 
search-and-rescue  operations. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


p 

pressure 

ADC 

analog-to-digital  convertor 

AWGN 

additive  white  gaussian  noise 

BER 

bit  error  rate 

DAC 

digital-to-analog  converter 

FDTD 

finite-difference  time-domain 

GNU 

GNU’s  Not  Unix 

GPS 

global  positioning  system 

HF 

high  frequency 

FOS 

line-of-sight 

MF 

matched  filter 

MUF 

MATE  AB  -USRP-Fibrary 

NLOS 

lon-line-of-lite 

QPSK 

ROS 

quadrature  phase  shift  keying 
robot  operating  system 

SDR 

software-defined  radio 

SINCGARS 

Single  Channel  Ground  and  Airborne  Radio  System 

SEAM 

simultaneous  localization  and  mapping 

SSD 

solid-state  drive 

UHF 

ultra  high  frequency 

USRP 

Universal  Software  Radio 

VHF 

very  high  frequency 

36 


Intentionally  lelt  blank. 


37 


1  DEFENSE  TECHNICAL 
(PDF)  INFORMATION  CTR 

DTIC  OCA 

2  DIRECTOR 

(PDF)  US  ARMY  RESEARCH  LAB 
RDRL  CIO  LL 

IMAL  HRA  MAIL  &  RECORDS  MGMT 

1  GOVT  PRINTG  OFC 

(PDF)  A  MALHOTRA 

8  US  ARMY  RESEARCH  LABORATORY 
(PDF)  RDRL-CIN-T 
F  DAGEFU 
G  VERMA 
C  RAO 
P  YU 
J  FINK 
B  SADLER 
B  RIVERA 
AKOTT 


38 


